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Abstract. We have performed detailed numerical calculations 
of the non-conservative evolution of close binary systems with 
low-mass (1.0 — 2.0M Q ) donor stars and a 1.3 M Q accreting 
neutron star. Rather than using analytical expressions for sim- 
ple polytropes, we calculated the thermal response of the donor 
star to mass loss, in order to determine the stability and follow 
the evolution of the mass transfer. Tidal spin-orbit interactions 
and Reimers wind mass-loss were also taken into account. 
We have re-calculated the correlation between orbital period 
and white dwarf mass in wide binary radio pulsar systems. Fur- 
thermore, we find an anti-correlation between orbital period 
and neutron star mass under the assumption of the "isotropic 
re-emission" model and compare this result with observations. 
We conclude that the accretion efficiency of neutron stars is 
rather low and that they eject a substantial fraction of the trans- 
ferred material even when accreting at a sub-Eddington level. 

The mass-transfer rate is a strongly increasing function of 
initial orbital period and donor star mass. For relatively close 
systems with light donors (P 01 ±, < 10 days and M% < 1.3 M®) 
the mass-transfer rate is sub-Eddington, whereas it can be 
highly super- Eddington by a factor of ~ 10 4 for wide systems 
with relatively heavy donor stars (1.6 ~ 2.0 M ) as a result of 
their deep convective envelopes. We briefly discuss the evolu- 
tion of X-ray binaries with donor stars in excess of 2 M@. 

Based on our calculations we present evidence that PSR 
J 1603-7202 evolved through a phase with unstable mass trans- 
fer from a relatively heavy donor star and therefore is likely to 
host a CO white dwarf companion. 
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1. Introduction 

Millisecond pulsars are characterized by short rotational pe- 
riods (P sp in < 30 ms) and relatively weak surface magnetic 
fields (B < 10 10 G) and are often found in binaries with a 
white dwarf companion. They are old neutron stars which have 
been recycled in a close binary via accretion of mass and angu- 
lar momentum from a donor star. The general scenario of this 
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process is fairly well understood qualitatively (cf. review by 
Bhattacharya & van den Heuvel 1991), but there remain many 
details which are still uncertain and difficult to analyze quanti- 
tatively. It is our aim to highlight these problems in a series of 
papers and try to answer them using detailed numerical calcu- 
lations with refined stellar evolution and binary interactions. 

There are now more than 30 binary millisecond pulsars 
known in the Galactic disk. They can be roughly divided into 
three observational classes (Tauris 1996). Class A contains the 
wide-orbit (P or b > 20 days) binary millisecond pulsars (BM- 
SPs) with low-mass helium white dwarf companions (1%d < 
0.45 M Q ), whereas the close-orbit BMSPs (P orb < 15 days) 
consist of systems with either low-mass helium white dwarf 
companions (class B) or systems with relatively heavy CO 
white dwarf companions (class C). The latter class evolved 
through a phase with significant loss of angular momentum 
(e.g. common envelope evolution) and descends from systems 
with a heavy donor star: 2 < M2/M0 < 6. The single millisec- 
ond pulsars are believed to originate from tight class B systems 
where the companion has been destroyed or evaporated - either 
from X-ray irradiation when the neutron star was accreting, or 
in the form of a pulsar radiation/wind of relativistic particles 
(e.g. Podsiadlowski 1991; Tavani 1992). 

The evolution of a binary initially consisting of a neutron 
star and a main-sequence companion depends on the mass of 
the companion (donor) star and the initial orbital period of the 
system. If the donor star is heavy compared to the neutron star 
then the mass transfer is likely to result in a common envelope 
(CE) evolution (Paczynski 1976; Webbink 1984; Iben & Livio 
1993) where the neutron star spirals in through the envelope 
of the donor in a very short timescale of less than 10 4 yr. The 
observational paucity of Roche-lobe filling companions more 
massive than ~ 2 Mq has been attributed to their inability to 
transfer mass in a stable mode such that the system becomes 
a persistent long-lived X-ray source (van den Heuvel 1975; 
Kalogera & Webbink 1996). For lighter donor stars (< 2 M ) 
the system evolves into a low-mass X-ray binary (LMXB) 
which evolves on a much longer timescale of 10 7 — 10 9 yr. It has 
been shown by Pylyser & Savonije (1988,1989) that an orbital 
bifurcation period (Pbif) separates the formation of converg- 
ing systems (which evolve with decreasing orbital periods until 
the mass-losing component becomes degenerate and an ultra- 
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compact binary is formed) from the diverging systems (which 
finally evolve with increasing orbital periods until the mass los- 
ing star has lost its envelope and a wide detached binary is 
formed). It is the LMXBs with P orb > Phu (— 2 days) which 
are the subject of this paper - the progenitors of the wide-orbit 
class A BMSPs. 

In these systems the mass transfer is driven by the inte- 
rior thermonuclear evolution of the companion star since it 
evolves into a (sub)giant before loss of orbital angular mo- 
mentum dominates. In this case we get an LMXB with a giant 
donor. These systems have been studied by Webbink, Rappa- 
port & Savonije (1983), Taam (1983), Savonije (1987), Joss, 
Rappaport & Lewis (1987) and recently Rappaportet al. (1995) 
and Ergma, Sarna & Antipova (1998). For a donor star on the 
red giant branch (RGB) the growth in core-mass is directly re- 
lated to the luminosity, as this luminosity is entirely generated 
by hydrogen shell burning. As such a star, with a small compact 
core surrounded by en extended convective envelope, is forced 
to move up the Hayashi track its luminosity increases strongly 
with only a fairly modest decrease in temperature. Hence one 
also finds a relationship between the giant's radius and the mass 
of its degenerate helium core - almost entirely independent 
of the mass present in the hydrogen-rich envelope (Refsdal & 
Weigert 1971; Webbink, Rappaport & Savonije 1983). In the 
scenario under consideration, the extended envelope of the gi- 
ant is expected to fill its Roche-lobe until termination of the 
mass transfer. Since the Roche-lobe radius i?L only depends on 
the masses and separation between the two stars it is clear that 
the core-mass, from the moment the star begins Roche-lobe 
overflow, is uniquely correlated with the orbital period of the 
system. Thus also the final orbital period, P^ rh (2 ~ 10 3 days) 
is expected to be a function of the mass of the resulting white 
dwarf companion (Savonije 1987). It has also been argued that 
the core-mass determines the rate of mass transfer (Webbink, 
Rappaport & Savonije 1983). For a general overview of the 
evolution of LMXBs - see e.g. Verbunt (1990). 
In this study we also discuss the final post-accretion mass of 
the neutron star and confront it with observations and the con- 
sequences of the new theory for kaon condensation in the core 
of neutron stars which result in a very soft equation-of-state and 
a corresponding maximum neutron star mass of only ~ 1.5M Q 
(Brown & Bethe 1994). 

In Section 2 we briefly introduce the code, and in Sections 3 
and 4 we outline the orbital evolution and the stability criteria 
for mass transfer. We present the results of our LMXB calcula- 
tions in Section 5 and in Section 6 we discuss our results and 
compare with observations. Our conclusions are given in Sec- 
tion 7 and a summary table of our numerical calculations is 
presented in the Appendix. 

2. A brief introduction to the numerical computer code 

We have used an updated version of the numerical stellar evo- 
lution code of Eggleton. This code uses a self-adaptive, non- 
Lagrangian mesh-spacing which is a function of local pressure, 
temperature, Lagrangian mass and radius. It treats both con- 



vective and semi-convective mixing as a diffusion process and 
finds a simultaneous and implicit solution of both the stellar 
structure equations and the diffusion equations for the chemical 
composition. New improvements are the inclusion of pressure 
ionization and Coulomb interactions in the equation-of-state, 
and the incorporation of recent opacity tables, nuclear reaction 
rates and neutrino loss rates. The most important recent updates 
of this code are described in Pols et al. (1995; 1998) and some 
are summarized in Han, Podsiadlowski & Eggleton (1994). 
We performed such detailed numerical stellar evolution calcu- 
lations in our work since they should result in more realistic 
results compared to models based on complete, composite or 
condensed poly tropes. 

We have included a number of binary interactions in this 
code in order to carefully follow the details of the mass-transfer 
process in LMXBs. These interactions include losses of orbital 
angular momentum due to mass loss, magnetic braking, gravi- 
tational wave radiation and the effects of tidal interactions and 
irradiation of the donor star by hard photons from the accreting 
neutron star. 

3. The equations governing orbital evolution 

The orbital angular momentum for a circulai{|] binary is: 

M NS M 2 



M 



Qa 2 



(1) 



where a is the separation between the stellar components, A/ns 
and A/ 2 are the masses of the (accreting) neutron star and the 
companion (donor) star, respectively, M = A/ns + AT 2 and 
the orbital angular velocity, £1 = ^/GM/a 3 . Here G is the 
constant of gravity. A simple logarithmic differentiation of this 
equation yields the rate of change in orbital separation: 
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where the total change in orbital angular momentum is: 
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The first term on the right side of this equation gives the change 
in orbital angular momentum due to gravitational wave radia- 
tion (Landau & Lifshitz 1958): 



J, 



orb 



32 G 3 M NS M 2 M 



(4) 



where c is the speed of light in vacuum. The second term arises 
due to magnetic braking. This is is a combined effect of a 
magnetically coupled stellar wind and tidal spin-orbit coupling 
which tend to keep the donor star spinning synchronously with 
the orbital motion. Observations of low-mass dwarf stars with 
rotational periods in the range of 1 ~ 30 days (Skumanich 



We assume circular orbits throughout this paper - tidal effects 
acting on the near RLO giant star will circularize the orbit anyway on 
a short timescale of ~ 10 4 yr, cf. Verbunt & Phinney (1995). 
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1972) show that even a weak (solar-like) wind will slow down 
their rotation in the course of time due to interaction of the 
stellar wind with the magnetic field induced by the differen- 
tial rotation in the convective envelope. For a star in a close 
binary system, the rotational braking is compensated by tidal 
coupling so that orbital angular momentum is converted into 
spin angular momentum and the binary orbit shrinks. Based on 
this observed braking law correlation between rotational period 
and age, Verbunt & Zwaan (1981) estimated the braking torque 
and we find: 



jmb ~ -0.5 x 10- 28 f~ 2 IR % 



Jovh 



mb 



a 5 Mns M 2 



(5) 



where R 2 is the radius of the donor star, I is its moment of iner- 
tia and / m b is a constant of order unity (see also discussion by 
Rappaport, Verbunt & Joss 1983). In order to sustain a signif- 
icant surface magnetic field we required a minimum depth of 
2 conv > 0.065 Rq for the convective envelope (cf. Pylyser & 
Savonije 1988 and references therein). Since the magnetic field 
is believed to be anchored in the underlaying radiative layers 
of the star (Parker 1955), we also required a maximum depth 
of the convection zone: Z conv / R 2 < 0.80 in order for the pro- 
cess of magnetic braking to operate. These limits imply that 
magnetic braking operates in low-mass (M 2 1.5 Mq) stars 
which are not too evolved. 

The third term on the right side of eq. (3) describes possi- 
ble exchange of angular momentum between the orbit and the 
donor star due to its expansion or contraction. For both this 
term and the magnetic braking term we estimate whether or 
not the tidal torque is sufficiently strong to keep the donor star 
synchronized with the orbit. The tidal torque is determined by 
considering the effect of turbulent viscosity in the convective 
envelope of the donor on the equilibrium tide. When the donor 
star approaches its Roche-lobe tidal effects become strong and 
lead to synchronous rotation. The corresponding tidal energy 
dissipation rate was calculated and taken into account in the 
local energy balance of the star. The tidal dissipation term was 
distributed through the convective envelope according to the lo- 
cal mixing-length approximation for turbulent convection - see 
Appendix for further details. 

Since we present calculations here for systems with P or b > 
2 days, the most significant contribution to the overall change 
in orbital angular momentum is caused by loss of mass from 
the system. This effect is given by: 



Jorb ~ 1 + q M 2 



(6) 



Here q = M 2 /Mns is the mass ratio of the donor over 
the accretor and a, [3 and 8 are the fractions of mass lost 
from the donor in the form of a fast wind, the mass ejected 
from the vicinity of the neutron star and from a circumstel- 
lar coplanar toroid (with radius, a r = j 2 a), respectively 
- see van den Heuvel (1994a) and Soberman, Phinney & 
van den Heuvel (1997). The accretion efficiency of the neutron 
star is thus given by: e = 1 — a — [3 — S, or equivalently: 



where dM 2 < 0. Note, that these factors will also be functions 
of time as the binary evolve. Low-mass (1 — 2 M Q ) donor stars 
do not lose any significant amount material in the form of a 
direct wind - except for an irradiated donor in a very close 
binary system, or an extended giant donor evolving toward the 
tip of the RGB which loses a significant amount of material in 
a wind. For the latter type of donors we used Reimers' (1975) 
formula to calculate the wind mass-loss rate: 



L R 

M 2 wind = -4x 10- 13 7 7 RW-rj 1 Moyr" 1 

M 2 



(8) 



where the mass, radius and luminosity are in solar units and 
77kw is the mass-loss parameter. We assumed rjnw = 0.5 for 
our work - cf. Renzini (1981) and Sackmann, Boothroyd & 
Kraemer (1993) for discussions. The mass-loss mechanism in- 
volving a circumstellar toroid drains far too much orbital an- 
gular momentum from the LMXB and would be dynamical un- 
stable resulting in a runaway event and formation of a CE. Also 
the existence of binary radio pulsars with orbital periods of sev- 
eral hundred days exclude this scenario as being dominant. 
Hence, for most of the work in this paper we have a <C (3, 
and we shall assume 5 = 0, and for LMXBs with large mass- 
transfer rates the mode of mass transfer to consider is therefore 
the "isotropic re-emission" model. In this model all of the mat- 
ter flows over, in a conservative way, from the donor star to 
an accretion disk in the vicinity of the neutron star, and then a 
fraction, (3 of this material is ejected isotropically from the sys- 
tem with the specific orbital angular momentum of the neutron 
star. 

As mentioned above, since we present calculations here for 
systems with initial periods larger than 2 days, loss of angular 
momentum due to gravitational wave radiation and magnetic 
braking (requiring orbital synchronization) will in general not 
be very significant. 

3.1. The mass-transfer rate 

For every timestep in the evolution calculation of the donor 
star the mass-transfer rate is calculated from the boundary con- 
dition on the stellar mass: 



M 2 = -1 x 10 3 PS 



Rl 



M q yr- 1 



(9) 



where PS[x] = 0.5 [x + abs(x)] and Rl is the donor's Roche- 
radius given by (Eggleton 1983): 



0.49 q 



2/3 , 



0.6<? 2 / 3 + ln(l + g 1/3 ) 



(10) 



<9M NS 



-(l-a-0-8) dM 2 



(7) 



The orbital separation a follows from the orbital angular mo- 
mentum balance - see eqs (1) and (3). All these variables are 
included in a Henyey iteration scheme. The above expression 
for the mass-transfer rate is rather arbitrary, as is the precise 
amount of Roche-lobe overfill for a certain transfer rate; but 
the results are independent of the precise form as they are de- 
termined by the response of the stellar radius to mass loss. 
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4. Stability criteria for mass transfer 

The stability and nature of the mass transfer is very important 
in binary stellar evolution. It depends on the response of the 
mass-losing donor star and of the Roche-lobe - see Soberman, 
Phinney & van den Heuvel (1997) for a nice review. If the mass 
transfer proceeds on a short timescale (thermal or dynamical) 
the system is unlikely to be observed during this short phase, 
whereas if the mass transfer proceeds on a nuclear timescale 
it is still able to sustain a high enough accretion rate onto the 
neutron star for the system to be observable as an LMXB for 
an appreciable interval of time. 

When the donor evolves to fill its Roche-lobe (or alterna- 
tively, the binary shrinks sufficiently as a result of orbital an- 
gular momentum losses) the unbalanced pressure at the first 
Lagrangian point will initiate mass transfer (Roche-lobe over- 
flow, RLO) onto the neutron star. When the donor star is per- 
turbed by removal of some mass, it falls out of hydrostatic and 
thermal equilibrium. In the process of re-establishing equilib- 
rium, the star will either grow or shrink - first on a dynami- 
cal (sound crossing), and then on a slower thermal (heat dif- 
fusion, or Kelvin-helmholtz) timescale. Also the Roche-lobe 
changes in response to the mass transfer/loss. As long as the 
donor star's Roche-lobe continues to enclose the star the mass 
transfer is stable. Otherwise it is unstable and proceeds on a 
dynamical timescale. Hence the question of stability is deter- 
mined by a comparison of the exponents in power-law fits of 
radius to mass, R oc M^, for the donor star and the Roche- 
lobe, respectively: 



Cdonor — 



<91ni? 2 
dlnM 2 



A Cl- 



d\nR L 
d\nM 2 



(11) 



where R 2 and M 2 refer to the mass losing donor star. Given 
R2 = Rl (the condition at the onset of RLO) the initial stabil- 
ity criterion becomes: 



CL < Cdo 



where (donor is the adiabatic or thermal (or somewhere in be- 
tween) response of the donor star to mass loss. Note, that the 
stability might change during the mass-transfer phase so that 
initially stable systems become unstable, or vice versa, later in 
the evolution. The radius of the donor is a function of time and 
mass and thus: 



#2 = 



dR 2 



dt 



8Rs. 



+ R2 Cdonor 



M2 

M 2 



dt 



M 2 

+ Rl Cl ~tt~ 

M 2 M 2 



(13) 



(14) 



The second terms follow from eq. (1 1); the first term of eq. (13) 
is due to expansion of the donor star as a result of nuclear burn- 
ing (e.g. shell hydrogen burning on the RGB) and the first term 
of eq. (14) represents changes in i?L which are not caused by 
mass transfer such as orbital decay due to gravitational wave 
radiation and tidal spin-orbit coupling. Tidal coupling tries to 



synchronize the orbit whenever the rotation of the donor is per- 
turbed (e.g. as a result of magnetic braking or an increase of 
the moment of inertia while the donor expands). The mass-loss 
rate of the donor can be found as a self-consistent solution to 
eqs (13) and (14) assuming R 2 = Rl for stable mass transfer. 

4.1. The Roche-radius exponent, Cl 

For binaries with orbital periods larger than a few days it is a 
good approximation that J gwr , J m b <C J m i and a <C f3 dur- 
ing the RLO mass-transfer phase. Assuming J gwr = J m b = 
and a = S = we can therefore use the analytical expression 
obtained by Tauris (1996) for an integration of eq. (2) to cal- 
culate the change in orbital separation during the LMXB phase 
(assuming a constant (3): 



3/3-5 

a_ = / go + i-/3 

a \q(l-0) + l 



1 



1 



(15) 



where the subscript '0' denotes initial values. Here we have 
added an extra factor, F\ s : 



Tig = cxp 2 



(dJ) ls 

Joib 



(16) 



to account for the tidal spin-orbit coupling since J\ s ^ 0. One 
aim of this study is to evaluate the deviation of T\ s from unity. 

If we combine eqs (7), (10) and (15), assuming Y\ s = 1, we 
obtain analytically: 



Cl 



(91ni? L fdlna <91n(_R L /a)\ dlnq 



dlnM 2 \dlnq 



dlnq 
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= [1 + (1 - (3) q] j> + (5 - 3/3) q 
where 

4 q 2/5 + i/3<r 1/3 (i + <7 1/3 ;r 1 



l + q 



0.6 + g- 2 / 3 ln(l + gV3) 



(17) 



(18) 



In the limiting case where q — > (when the accretor is much 
heavier than the donor star) we find: 



lim Cl 

q^O 



-5/3 



(19) 



The behavior of Cl (q, (3) for LMXBs is shown in Fig. 1. We 
note that Cl does not depend strongly on (3. This figure is quite 
useful to get an idea of the stability of a given mass trans- 
fer when comparing with C for the donor star. We see that in 
general the Roche-lobe, _Rl increases (Cl < 0) when mate- 
rial is transferred from a light donor to a heavier NS (q < 1) 
and correspondingly i?L decreases (Cl > 0) when material 
is transferred from a heavier donor to a lighter NS (q > 1). 
This behavior is easily understood from the bottom panel of 
the same figure where we have plotted — d\n(a) /d\n(q) as a 
function of q. The sign of this quantity is important since it 
tells whether the orbit expands or contracts in response to mass 
transfer (note dq < 0). We notice that the orbit always expands 
when q < 1 and it always decreases when q > 1.28, whereas 
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Fig. 1. Top panel: the Roche-radius exponent (i?L oc M%) for 
LMXBs as a function of q and j3. The different curves corre- 
spond to different constant values of (3 in steps of 0.1. Tidal 
effects were not taken into account (T\ s = 1). A cross is shown 
to highlight the case of q = 1 or £l = 0. In the bottom panel 
we have plotted — d ln(a)/d hi(q) as a function of q. The evo- 
lution during the mass-transfer phase follows these curves from 
right to left (though (3 need not be constant) since M 2 and q are 
decreasing with time. See text for further explanation. 

for 1 < q ^ 1.28 it can still expand if f3 > 0. There is a point 
at q = 3/2 where d ln(a) / d ln(q) — 2/5 is independent of f3. 
It should be mentioned that if j3 > then, in some cases, it 
is actually possible to decrease the separation, a between two 
stellar components while increasing P or \, at the same time! 

5. Results 

We have evolved a total of a few hundred different LMXB 
systems. 121 of these are listed in Table Al in the Appendix. 
We chose donor star masses of 1.0 < M%/Mq < 2.0 and 
initial orbital periods of 2.0 < P Q z r AMS /days < 800. We 
also evolved donors with different chemical compositions and 
mixing-length parameters. In all cases we assumed an initial 
neutron star mass of 1.3 Mq. 

In Fig. 2 we show the evolution of four LMXBs. As a function 



of donor star mass (M2) or its age since the ZAMS, we have 
plotted the orbital period (P or b)> the mass-loss rate of the donor 
as well as the mass-loss rate from the system (A/2 and M), the 
radius exponent (Q of the donor and its Roche-lobe and finally 
the depth of the donor's convection zone (Z com j R2). Note, 
that we have zoomed in on the age interval which corresponds 
to the mass-transfer phase. As an example, we have chosen two 
different initial donor masses (1.0 Mq and 1.6 Mq) — each with 
two different initial orbital periods (3.0 and 60.0 days) of the 
neutron star (NS) and its ZAMS companion. The evolutionary 
tracks of the donor stars are plotted in the HR-diagram in Fig. 3. 
We will now discuss the evolution of each of these systems in 
more detail. 



and P ZAMS = 3.0 



5.7. Fig. 2a 

In Fig. 2a we adopted Af| = 1.0 Mq 
days. In this case the time it takes for the donor to become 
a (sub)giant and fill its Roche-lobe, to initiate mass transfer, 
is 11.89 Gyr. Before the donor fills its Roche-lobe the expan- 
sion due to shell hydrogen burning causes its moment of inertia 
to increase which tends to slow down the rotation of the star. 
However, the tidal torques act to establish synchronization by 
transferring angular momentum to the donor star at the expense 
of orbital angular momentum. Hence at the onset of the mass 
transfer (A) the orbital period has decreased from the initial 
P^ams _ 3 days to P^° = 1.0 days and the radius is now 
R2 = Rh = 2.0 Rq. 

We notice that the mass-loss rate of the donor star remains 
sub-Eddington (\M 2 \ < M Edd rj 1.5 x 10" 8 Af Q yr- 1 for 
hydrogen-rich matter) during the entire mass transfer]^. Thus 
we expect all the transferred material to be accreted onto the 
neutron star, if disk instabilities and propeller effects can be ne- 
glected (see Sections 5.7 and 6.4). Therefore we have no mass 
loss from the system in this case - i.e. M = 0. The duration 
of the mass-transfer phase for this system is quite long: ~1.0 
Gyr (11.89 -> 12.91 Gyr). At age, t ~ 12.65 Gyr (P orb = 5.1 
days; M 2 = 0.317 Mq) the donor star detaches slightly from 
its Roche-lobe (d) and the mass transfer ceases temporarily for 
~ 25 Myr - see next subsection for an explanation. 
The Roche-radius exponent calculated from eq. (17) is plot- 
ted as a dotted line as a function of M2 in the upper right 
panel. However, our numerical calculations (full line) show that 
tidal effects are significant and increase ( by ~ 0.5-0.8 until 
Mi rs 0.54 Mq (p). At this point the magnetic braking is as- 
sumed to switch off, since Z conv /i?2 > 0.80. Note that during 
the mass transfer phase £ ss £ L and, as long as the mass trans- 
fer is not unstable on a dynamical timescale, we typically have 
in our code: 1 x 10" 4 < ln(R 2 /Ri J ) < 7 x 10~ 3 and hence 



2 Strictly speaking MEdd = AfEdd(7?Ns) is slightly reduced dur- 
ing the accretion phase since the radius of the neutron star decreases 
with increasing mass (e.g. for an ideal n-gas polytrope: i?NS oc 

X /3 

M NS ). However, this only amounts to a correction of less than 20% 
for various equations-of-state, and thus we have not taken this effect 
into account. 
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M 2 (M ) Age (yr) M 2 (M ) 



Fig.2a-b. Numerical calculations of the mass-transfer process in two LMXB systems with -P^b^ 13 = ^-0 days, X = 0.70, Z = 
0.02 and a = 2.0. The mass of the donor stars are (a): 1.0 M Q , and (b): 1.6 M Q - top and bottom panels, respectively. See text 
for further details. 
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practically C = Cl- 

The final outcome for this system is a BMSP with an orbital 
period of P^ rb = 9.98 days and a He white dwarf (WD) with 
a mass of M WD = 0.245 M© (B). The final mass of the NS 
is Mns = 2.06 Mq, since we assumed all the material was 
accreted onto the NS given \M 2 \ < A^Edd during the entire 
X-ray phase. However, in Section 6 we will discuss this as- 
sumption and the important question of disk instabilities and 
the propeller mechanism in more detail. 

5.2. Fig. 2b 

In Fig. 2b we adopted M\ = 1.6 M© and P z r £ MS = 3.0 

days. The RLO is initiated at an age of t = 2.256 Gyr when 
pRLO = 2 days and R2 = 3 8 Rq (A) In this 

case the mass- 
transfer rate is super- Eddington (\M 2 \ > MEdd, cf. dashed 
line) at the beginning of the mass-transfer phase. In our adopted 
model of "isotropic re-emission" we assume all material in ex- 
cess of the Eddington accretion limit to be ejected from the 
system, while carrying with it the specific orbital angular mo- 
mentum of the neutron star. Hence \M\ = \M 2 \ — M E dd- Ini- 
tially \M 2 \ ~ 10 2 M E dd at the onset of the RLO and then \M 2 \ 
decreases from 10 to 1 Af E dd at M 2 ~ 0.7M© (e) before it 
becomes sub-Eddington for the rest of the mass-transfer phase. 
Mass loss from the system as a result of a Reimers wind in 
the red giant stage prior to RLO (A) is seen to be less than 
10~ n M© yr _1 . By comparing the different panels for the evo- 
lution, we notice that the initial super-Eddington mass transfer 
phase ( A — e) lasts for 22 Myr. In this interval the companion 
mass decreases from 1.6 Mq to 0.72 Mq. Then the system en- 
ters a phase (e—d) of sub-Eddington mass transfer at P or b=5.3 1 
days which lasts for 41 Myr. When M 2 = 0.458 M©, and 
-Porb=13.6 days, the system detaches and the X-ray source is 
extinguished for about 40 Myr (d), cf. gray-shaded area. The 
temporary detachment is caused by a transient contraction of 
the donor star when its hydrogen shell source moves into the 
hydrogen rich layers left behind by the contracting convective 
core during the early main sequence stage. At the same time 
the convective envelope has penetrated inwards to its deepest 
point, i.e. almost, but not quite, to the H-shell source. The effect 
of a transient contraction of single low-mass stars evolving up 
the RGB, as a result of a sudden discontinuity in the chemical 
composition, has been known for many years (Thomas 1967; 
Kippenhahn & Weigert 1990) but has hitherto escaped atten- 
tion in binary evolution. After the transient contraction the star 
re-expands enough to fill its Roche-lobe again and further as- 
cends the giant branch. The corresponding final phase of mass 
transfer (d-B) is sub-Eddington (|M 2 | « 0.2 M E dd) and lasts 
for 60 Myr. The end product of this binary is a recycled pulsar 
and a He- WD companion with an orbital period of 41 .8 days. In 
this case we obtain M WD = 0.291 M and M NS = 2.05 Mq. 
The total duration of the mass-transfer phase during which the 
system is an active X-ray source is tx = 123 Myr (excluding 
the quiescence phase of 40 Myr) which is substantially shorter 
compared to the case discussed above (Fig. 2a). 
The reason for the relatively wide final orbit of this system, 
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Fig. 3. Calculated evolutionary tracks in the Hertzsprung- 
Russel diagram of the four donor stars used in Fig. 2. The 
dashed lines represent the evolution of a single star with mass 
M 2 . The mass-transfer phase (RLO) is initiated at A and ceases 
at B. The symbols along the evolutionary tracks correspond to 
those indicated in Fig. 2. 



compared to the case discussed above with the 1.0 Mq donor, 
is caused by the super-Eddington mass transfer during which a 
total of 0.55 Mq is lost from the system. 
The numerical calculations of ( for this donor star (full line) fits 
very well with our simple analytical expression (dotted line) 
which indicates that the effects of the tidal spin-orbit interac- 
tions are not so significant in this case. 
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5.3. Fig. 2c 

In this figure we adopted M| = 1.0 Mq and P z r £ MS = 60.0 
days. The RLO is initiated (A) at an age of t ~ 12.645 
Gyr. At this stage the mass of the donor has decreased to 
y\,/RLO _ o.976 Mq as a result of the radiation-driven wind 
of the giant star. However, the orbital period has also decreased 
(■^OTb = ^8-1 days) and thus the shrinking of the orbit due 
to tidal spin-orbit coupling dominates over the widening of the 
orbit caused by the wind mass loss. 

The total interval of mass transfer is quite short, t% — 13.3 
Myr. The mass-transfer rate is super-Eddington during the en- 
tire evolution (IM2I ~ 1 — 6A/Edd) and therefore the NS 
only accretes very little material: AMns = (MEdd) At mt ~ 
0.20 Mq. The reason for the high mass-loss rate of the donor 
star is its deep convective envelope (see lower right panel). 
Since the initial configuration of this system is a very wide or- 
bit, the donor will be rather evolved on the RGB when it fills its 
Roche-lobe (R 2 = 29.3 R Q and = 58.1 days). Hence 

the donor swells up in response to mass loss (i.e. ( < 0) as 
a result of the super-adiabatic temperature gradient in its giant 
envelope. The radius exponent is well described by our ana- 
lytical formula in this case. The final outcome of this system 
is a wide-orbit (P Q f rb = 382 days) BMSP with a - 0.40 M@ 
He- WD companion. 

5.4. Fig. 2d 

Here we adopted M| = 1.6 M Q and -P z r £ MS = 60.0 days. At 
the onset of the RLO the donor mass is Af^ LO = 1.582 Mq. In 
this case we do not only have a giant donor with a deep convec- 
tive envelope. It is also (initially) heavier than the accreting NS. 
Both of these circumstances makes it difficult for the donor to 
retain itself inside its Roche-lobe once the mass transfer is ini- 
tiated. It is likely that such systems, with huge mass-transfer 
rates, evolve into a phase where matter piles up around the 
neutron star and presumably forms a growing, bloated cloud 
engulfing it. The system could avoid a spiral-in when it man- 
ages to evaporate the bulk of the transferred matter from the 
surface of the (hot) accretion cloud via the liberated accretion 
energy. This scenario would require the radius of the accre- 
tion cloud, r c \ to be larger than ~ 7?ns M-^dd) m order 
for the liberated accretion energy to eject the transfered mate- 
rial. However, if there is insufficient gas cooling r c \ could be 
smaller from an energetic point of view. At the same time r c \ 
must be smaller than the Roche-lobe radius of the neutron star 
(cf. eq. 10 with q — M-ns/Mz) during the entire evolution. In 
that case our simple isotropic re-emission model would approx- 
imately remain valid. Assuming this to be the case we find the 
mass-transfer rate is extremely high: \M 2 \ ~ 10 4 A^Edd and 
more than 0.5 Mq is lost from the donor (and the system) in 
only a few 10 3 yr. The system survives and the orbital period 
increases from 54.5 days to 1 1 1 days during this short phase. 
After this extremely short mass-transfer epoch, with an ultra- 
high mass-transfer rate, the donor star relaxes (d) and shrink 
inside its Roche-lobe for 2.5 Myr when M 2 = 0.98 Mq. The 



mass transfer is resumed again for 7.5 Myr at a more moderate 
super-Eddington rate (d — E). The final outcome is a binary 
pulsar with a 0.43 Mq He- WD companion and an orbital pe- 
riod of 608 days. Though the NS only accretes ~ 0.10 Mq as a 
result of the short integrated accretion phase it will probably be 
spun-up sufficiently to become a millisecond pulsar since mil- 
lisecond pulsars evidently are also formed in systems which 
evolve e.g. through a CE with similar (or even shorter) phases 
of accretion (van den Heuvel 1994b). 

The initial extreme evolution of this system causes an offset in 
Q until the more moderate mass-transfer phase (d — E) contin- 
ues at M% — 0.98 Mq. It should be noted that a system like this 
is very unlikely to be observed in the ultra-high mass-transfer 
state due to the very short interval (< 10 4 yr) of this phase. 

5.5. M2 > 2 Mq, runaway mass transfer and onset of a CE 

The latter example above illustrates very well the situation 
near the threshold for unstable mass transfer on a dynamical 
timescale and the onset of a CE evolution]]. If the donor star is 
heavier than 1.8 Mq a critical overflow is likely to occur since 
the orbit shrinks in response to mass transfer (q < 1.28, cf. 
Section 4). This is also the situation if P or b is large because the 
donor in that case develops a deep convective envelope which 
causes it to expand in response to mass loss and a runaway 
mass transfer sets in. When a runaway mass transfer sets in 
we were not able to prevent it from critically overflowing its 
Roche-lobe and our code breaks down. At this stage the neu- 
tron star is eventually embedded in a CE with its companion 
and it will spiral in toward the center of its companion as a re- 
sult of removal of orbital angular momentum by the drag force 
acting on it^|. The final result of the CE depends mainly on the 
orbital period and the mass of the giant's envelope. If there is 
enough orbital energy available (i.e. if P or b is large enough at 
the onset of the CE), then the entire envelope of the giant can 
be expelled as a result of the liberated orbital energy, which is 
converted into kinetic energy that provides an outward motion 
of the envelope decoupling it from its core. This leaves behind a 
tight binary with a heavy WD (the core of the giant) and a mod- 
erately recycled pulsar. There are five such systems observed in 
our Galaxy. They all have a CO- WD and P or b = 6^8 days. 
These are the so-called class C BMSPs. 

If there is not enough orbital energy available to expel the enve- 
lope, then the NS spirals in completely to the center of the giant 
and a Thorne-Zytkow object is formed. Such an object might 
evolve into a single millisecond pulsar (e.g. van den Heuvel 
1994a) or may collapse into a black hole (Chevalier 1996). 

3 We notice, that this very high mass-transfer rate might lead to 
hyper-critical accretion onto the neutron star and a possible collapse 
of the NS into a black hole if the equation-of-state is soft (cf. Cheva- 
lier 1993; Brown & Bethe 1994; Brown 1995). However, new results 
obtained by Chevalier (1996) including the centrifugal effects of a ro- 
tating infalling gas might change this conclusion. 

4 However, even binaries with donor stars of 2-6 Mq might survive 
the mass transfer avoiding a spiral-in phase in case the envelope of the 
donor is still radiative at the onset of the RLO. 
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5.6. The (P or b, Mwd) correlation 

We have derived new (P or b, Mwd) correlations based on the 
outcome of the 121 LMXB models calculated for this work. 
They are shown in Fig. 4a (the top panel). We considered mod- 
els with donor star masses (M 2 ) in the interval 1.0 — 2.0 M©, 
chemical compositions ranging from Population I (X=0.70, 
Z=0.02) to Population II (X=0.75, Z=0.001) and convective 
mixing-length parameters a = 1/ H p from 2-3 (here I is the 
mixing length and H p is the local pressure scaleheight). Fol- 
lowing Rappaport et al. (1995) we chose our standard model 
with M 2 = 1.0 M©, Population I composition and a = 2, cf. 
thick line in Fig. 4. The upper limit of M 2 is set by the require- 
ment that the mass transfer in the binary must be dynamically 
stable, and the lower limit by the requirement that the donor 
star must evolve off the main sequence within an interval of 
time given by: t ms < i Hu bbio - igai - t coo \. Here t Hu bbio ~ 15 
Gyr is the age of the Universe, i ga i ~ 1 Gyr is the minimum 
time between the Big Bang and formation of our Milky Way 
and t coo \ ~ 3 Gyr is a typical low value of WD companion 
cooling ages, following the mass-transfer phase, as observed in 
BMSPs (Hansen & Phinney 1998). We thus find M 2 ~ 1.0 M 
as a conservative lower limit. 

The first thing to notice, is that the correlation is more or 
less independent of the initial donor star mass (M 2 ) - only for 
M 2 <; 2.0M© (where the mass transfer becomes dynamically 
unstable anyway for P* lb <; 4.2 days) we see a slight devia- 
tion. This result is expected if M 2 coro (and therefore R 2 and 
Porb) is independent of M 2 . We have performed a check on 
this statement using our calculations for an evolved donor star 
on the RGB. As an example, in Table 1 we have written L, T err 
and M 2 coro as a function of M 2 when it has evolved to a ra- 
dius of 50.0 Rq. In addition we have written the mass of the 
donor's envelope at the moment R 2 = 50.0 Rq. We conclude, 

Table 1. Stellar parameters for a star with R 2 = 50.0 R Q - see 
text. 



M2/M0 


1.0** 


1.6** 


1.0* 


1.6* 


logi/L© 


2.566 


2.624 


2.644 


2.723 


logTrff 


3.554 


3.569 


3.573 


3.593 


M 2c ore/M 


0.336 


0.345 


0.342 


0.354 


M 2cn v/M 


0.215 


0.514 


0.615 


1.217 



* Single star (X=0.70, Z=0.02 and a=2.0). 
** Binary donor (P z r £ MS = 60.0 days and M NS 



1.3 M Q ) 



for a given chemical composition and mixing-length parameter, 
M 2 corc is practically independent of M 2 (to within a few per 
cent) and that mass loss from the envelope via RLO has simi- 
lar little influence on the (R 2 , M 2 coro ) correlation as well. For 
other choices of R 2 the differences were found to be smaller. In 
Fig. 5 we have shown P or b (which resembles R 2 ) as a function 
of M 2coro . 

Much more important is the theoretical uncertainty in the value 



of the convective mixing-length parameter and most important 
is the initial chemical composition of the donor star. We have 
estimated a (-P rb> -Mwd) correlation from an overall best fit 
to all the models considered in Table Al and obtain (P or b in 
days): 



Mi 



WD 



1/a 



+ C 



(20) 



Mrs, V b 

where, depending on the chemical composition of the donor, 

( 4.50 1.2 x 10 5 0.120 Pop.I 
(a,b,c) = } 4.75 1.1 x 10 5 0.115 Pop.I+II (21) 

[ 5.00 1.0 x 10 5 0.110 Pop.II 

This formula is valid for binaries with: 0.18 ^ M^,° D /M Q <; 
0.45. The uncertainty in the initial chemical abundances of the 
donor star results in a spread of a factor ~1.4 about the me- 
dian (Pop.I+II) value of P or b at any given value of Mwd- The 
spread in the (P or b,MwD) correlation arises solely from the 
spread in the (R 2 , M 2colc ) correlation as a result of the dif- 
ferent chemical abundances, and/or a, of the giant donor star 
ascending the RGB. 

If we compare our calculations with the work of Rappa- 
port et al. (1995) we find that our best fit results in signifi- 
cantly lower values of P or b for a given mass of the WD in 
the interval 0.18 ^ M W d/M q ^ 0.35. It is also notable 
that these authors find a maximum spread in P or b of a factor 
-2.4 at fixed M WD . For 0.35 £ M WD /M© £ 0.45 their 
results agree with our calculations to within 20%. A fit to 
our eq. (20) with the results of Rappaport et al. (1995) yields: 
a = 5.75, b = 8.42 x 10 4 and c = (to an accuracy within 
1 % for 0.18 <> M WD /M© & 0.45) for their donor models 
with population I chemical composition and a = 2.0. For their 
Pop. II donors we obtain b = 3.91 x 10 4 and same values for 
a and c as above. We also obtain somewhat lower values of 
Porb, for a given mass of the WD, compared with the results of 
Ergma, Sarna & Antipova (1998). 

The discrepancy between the results of the above men- 
tioned papers and our work is a result of different input physics 
for the stellar evolution (cf. Section 2). Ergma, Sarna & An- 
tipova (1998) uses models based on Paczynzki's code, and Rap- 
paport et al. (1995) used an older version of Eggleton's code 
than the one used for this work. In our calculations we have 
also included the effects of tidal dissipation. However, these 
effects can not account for the discrepancy since in this paper 
we only considered binaries with P orb > 2 days and thus the 
effects of the tidal forces are relatively small (the contribution 
to the stellar luminosity from dissipation of tidal energy is only 
itidai ^ 0.05L nuc for P orb = 2 days). 

In analogy with Rappaport et al. (1995) and Ergma, Sarna 
& Antipova (1998) we find that, for a given value of M WD , 
P or b is decreasing with increasing a, and P or b is increasing 
with increasing metallicity. We find M^ 001 ~ M^ 02 + 
0.03M© which gives a stronger dependency on metallicity, by 
a factor —2, compared to the work of Ergma, Sarna & Antipova 
(1998). 
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Fig. 4a-b. The (P rb, Mwd) correlation (top panel) and the (P rb, A^ns) anti-correlation (bottom panel) calculated for different 
donor star masses, M2 chemical compositions (Pop.l: X=0.70; Z=0.02 and Pop. II: X=0.75; Z=0.001) and mixing-length param- 
eters, a as indicated in the top panel. The gray line shows the correlation obtained by Rappaport et al. (1995) for M2 = 1.0 Mq, 
Pop.l chemical abundances and a = 2.0. The post-accretion Mns curves (bottom) assume no mass loss from accretion disk 
instabilities of propeller effects - see Sections 5.7 and 6.4. 
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Fig. 5. Evolutionary tracks of the four LMXBs in Figs 2 and 
3, showing the binary orbital period changing as a function of 
the mass of the core of the donor star. At the termination of 
the mass-transfer process M 2cor c w Mwd — 0.01 M Q and 
the end-points of the evolutionary tracks are located near the 
curve representing the (P or b, Mwd) correlation. The initial or- 
bital periods were P z r £ MS = 3.0 days and P Z £ MS = 60.0 days 
for the two bottom and top tracks, respectively. Furthermore we 
used Population I chemical abundances and a = 2. 

It should be noticed, that the (P or b, Mwd) correlation is in- 
dependent of (3 (the fraction of transferred material lost from 
the system), the mode of mass loss and degree of magnetic 
braking since, as demonstrated above, the relationship between 
i? 2 and M 2 CO ro of the giant donors remains unaffected by the 
exterior stellar conditions governing the process of mass trans- 
fer - see also Alberts et al. (1996). But for the individual binary, 
P or b and Mwd do depend on (3 and they increase with increas- 
ing values of j3 (see e.g. the bottom panel of Fig. 1 for q < 1 
which always applies near the end of the mass-transfer phase). 

As mentioned in our examples earlier in this section, there 
is a competition between the wind mass loss and the tidal spin- 
orbit interactions for determining the orbital evolution prior to 
the RLO-phase. This is demonstrated in Fig. 6 where we have 
shown the changes in P or b and M 2 , from the ZAMS stage to 
the onset of the RLO, as a function of the initial ZAMS orbital 
period. It is seen that only for binaries with P z r £ MS > 100 
days will the wind mass-loss be efficient enough to widen the 
orbit. For shorter periods the effects of the spin-orbit interac- 
tions dominate (caused by expansion of the donor) and loss of 
orbital angular momentum causes the orbit to shrink. 

5.7. The ( P or b, Mns) anti-correlation 

We now investigate the interesting relationship between the fi- 
nal mass of the NS and the final orbital period. In Fig. 4b (the 
bottom panel) we have plotted P or b as a function of the poten- 
tial maximum mass of the recycled pulsar, M^ s . This value is 




n ZAMS / i \ 

Fig. 6. The changes of donor mass, M 2 (full lines) and orbital 
period, P or b (dashed lines), due to wind mass loss and tidal 
spin-orbit interactions, from the ZAMS until the onset of the 
RLO as a function of the initial orbital period. Plots are shown 
for two different values of the Reimers' mass-loss parameter, 
77rw- The binary is assumed to be circular. See text for further 
discussion. 



the final mass of the NS if mass loss resulting from instabil- 
ities in the accretion process are neglected. Another (smaller) 
effect which has also been ignored is the mass deficit of the ac- 
creted material as it falls deep down the gravitational potential 
well of the NS. The gravitational mass of a NS (as measured 
from a distant observer by its gravitational effects) contains not 
only the rest mass of the baryons, but also the mass equiva- 
lent of the negative binding energy, AMdcf = Pbind/c 2 < 0. 
Depending on the equation-of-state AM dol ^10% of the grav- 
itational mass (Shapiro & Teukolsky 1983). This is hence also 
the efficiency of radiative emission in units of available rest- 
mass energy incident on the NS. Thus we can express the ac- 
tual post-accretion gravitational mass of a recycled pulsar by 
(<9m 2 < 0): 



M NS = M^ s + 



JM 2 



Mwd 



(1 - /?') dm 2 - AM dp 



fcdcf(22) 



Here j3' = max ^(|M 2 | - M Edd )/|M 2 |, 0J is the fraction of 
material lost in a relativistic jet as a result of super-Eddington 
mass transfer; AM dp = AM disk + AM prop is the sum of mat- 
ter lost from the accretion disk (as a result of viscous instabil- 
ities or wind corona) and matter being ejected near the pulsar 
magnetosphere as a result of the centrifugal propeller effect, 

and finally fc dof = ( M NS ^AM d f ) ~ °- 90 is a factor that ex- 
presses the ratio of gravitational mass to rest mass of the ma- 
terial accreted onto the NS. M^ s used in Fig. 4b is given by 
the expression above assuming AM dp = and AM dcf = 
(fc def = 1). 
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We see that the (P or b, M^ s ) anti-correlation is more or less 
independent of the chemical composition and a of the donor 
star, whereas it depends strongly on M 2 for P orb & 50 days. 
This anti-correlation between P or b and M^ s is quite easy to 
understand: binaries with large initial orbital periods will have 
giant donor stars with deep convective envelopes at the onset 
of the mass transfer; hence the mass-transfer rate will be super- 
Eddington and subsequently a large fraction of the transferred 
material will be lost from the system. Therefore BMSPs with 
large values of P or b are expected to have relatively light NS 
- cf. Sections 5.3 (Fig. 2c) and 5.4 (Fig. 2d). Similarly, bina- 
ries with small values of Po r t MS will result in BMSPs with 
relatively small P or b and large values of M^ s , since M2 will 
be sub-Eddington and thus the NS has the potential to accrete 
all of the transferred material - cf. Sections 5.1 (Fig. 2a) and 
5.2 (Fig. 2b). Therefore, if disk instabilities, wind corona and 
propeller effects were unimportant we would expect to find 
an (P or b,M N s) anti-correlation among the observed BMSPs. 
However, below (Section 6.2) we demonstrate that mass ejec- 
tion arising from these effects is indeed important and thus it is 
very doubtful whether an (P or b, Mns) anti-correlation will be 
found from future observations. 



6. Discussion and comparison with observations 

6.1. Comparison with condensed polytrope donor models 

Hjellming & Webbink (1987) studied the adiabatic properties 
of three simple families of polytropes by integrating the non- 
linear Lane-Emden equation in Lagrangian coordinates. The 
condensed polytropes, consisting of n = 3/2,7 = 5/3 (con- 
vective) envelopes with He-core point masses, are suitable for 
red giant stars. In is not trivial to directly compare our calcu- 
lations with e.g. the stability analysis of Soberman, Phinney & 
van den Heuvel (1997) and Kalogera & Webbink (1996) since 
the donor does not restore (thermal) equilibrium after initia- 
tion of an unstable mass-transfer process. But it is important to 
point out, that systems which initiate RLO with thermally un- 
stable mass transfer could, in some cases, survive this tempo- 
rary phase - even if | M2 1 exceeds the Eddington accretion limit 
by as much as a factor ~10 4 (see Fig. 2d). Similarly, systems 
which begin mass transfer on a thermal timescale may in some 
cases (if M2 is large compared to Mns) eventually become dy- 
namically unstable. These results were also found by Hjellming 
(1989) and Kalogera & Webbink (1996), and we refer to these 
papers for a more detailed discussion on the fate of thermally 
unstable systems. Therefore it is not always easy to predict the 
final outcome of an LMXB system given its initial parameters 
- especially since the onset criteria of a CE phase is rather 
uncertain. Nevertheless, we can conclude that LMXBs with 
M2 < 1.8 Mq will always survive the mass transfer phase. 
Systems with donor stars M2 > 2 M Q only survive if P^ MS 
is within a certain interval. 

Soberman, Phinney & van den Heuvel (1997) also used 
the polytropic models of Hjellming & Webbink (1987) to fol- 
low the mass transfer in binaries. The global outcome of such 



calculations are reasonably good. However, the weakness of 
the polytropic models is that whereas they yield the radius- 
exponent at the onset of the mass transfer, and the approxi- 
mated stellar structure at that given moment, they do not trace 
the response of the donor very well during the mass-transfer 
phase. The structural changes of the donor star (e.g. the out- 
ward moving H-shell and the inward moving convection zone 
giving rise to the transient detachment of the donor from its 
Roche-lobe) can only be followed in detail by a Henyey-type 
iteration scheme for a full stellar evolutionary model. 

6.2. The observed (Porb, Mwd) correlation 

The companion mass, Mwd of an observed binary pulsar is 
constrained from its Keplerian mass function which is obtained 
from the observables P or b and a p sin i: 



f (M NS , M WD ) 



(M WD sini) 3 4.7T 2 (a p sini) 3 



(Mns + Mwd) 2 



G 



Pi 



(23) 



Here i is the inclination angle (between the orbital angular mo- 
mentum vector and the line-of-sight toward the observer) and 
a p = a (Mwd/(Mns + Mwd)) is the semi-major axis of the 
pulsar in a cm. reference frame. The probability of observing 
a binary system at an inclination angle i, less than some value 
io, is P(i < iq) = 1 — cos («o). 

As mentioned earlier, there are indeed problems with fitting 
the observed low-mass binary pulsars onto a theoretical core- 
mass period relation. The problem is particularly pronounced 
for the very wide-orbit BMSPs. Although the estimated masses 
of the companions are quite uncertain (because of the un- 
known orbital inclination angles and Mns) no clear observed 
(P or b, Mwd) correlation seems to be present - opposite to 
what is proposed by several authors (e.g. Phinney & Kulka- 
rni 1994, Lorimer et al. 1996 and Rappaport et al. 1995). In 
Table A2 in the Appendix we have listed all the observed 
galactic (NS+WD) binary pulsars and their relevant parame- 
ters. It was noticed by Tauris (1996) that the five BMSPs with 
P>rb > 100 days all seem to have an observed M^, which 
is lighter than expected from the theoretical correlation (at the 
^80 % confidence level on average). There does not seem to 
be any observational selection effects which can account for 
this discrepancy (Tauris 1996; 1998) - i.e. why we should pref- 
erentially observe systems with small inclination angles (sys- 
tematic small values of i, rather than a random distribution, 
would increase M^, for the given observed mass functions 
and thus the observations would match the theory). Evapora- 
tion of the companion star, from a wind of relativistic particles 
after the pulsar turns on, also seems unlikely since the evap- 
oration timescale (proportional to P^) becomes larger than 
^Hubble for such wide orbits. It is also worth mentioning that 
the orbital period change due to evaporation, or general mass- 
loss in the form of a stellar wind, is at most a factor of ^2, if 
one assumes the specific angular momentum of the lost matter 
is equal to that of the donor star. 

Beware that the (P or b, Mwd) correlation is not valid for BM- 
SPs with CO/O-Ne-Mg WD companions as these systems did 
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not evolve through a phase with stable mass transfer. The ex- 
ception here are the very wide orbit systems with P^ rb <; 800 
days. PSR B0820+02 might be an example of such a system. 
From Table Al (Appendix) it is seen that we expect a maxi- 
mum orbital period of ~ 1400 days for the NS+WD binaries. 
Larger periods are, of course, possible but the binaries are then 
too wide for the neutron star to be recycled via accretion of 
matter. 

It should also be mentioned that the recycling process is 
expected to align the spin axis of the neutron star with the or- 
bital angular momentum vector as a result of > 10 7 yr of stable 
disk accretion. Hence we expect (Tauris 1998) the orbital in- 
clination angle, i to be equivalent to (on average) the magnetic 
inclination angle, a mag defined as the angle between the pulsar 
spin axis and the center of the pulsar beam (viz. line-of-sight to 
observer). 

6.3. PSRJ2019+2425 

PSR J20 19+2425 is a BMSP with P orb = 76.5 days and a mass 
function / = 0.0107M© (Nice, Taylor & Fruchter 1993). In 
a recent paper (Tauris 1998) it was demonstrated that for this 
pulsar Mns — 1.20 M Q , if the (P or b,MwD) correlation ob- 
tained by Rappaport et al. (1995) was taken at face value. This 
value of Mns is significantly lower than that of any other es- 
timated pulsar mass (Thorsett & Chakrabarty 1999). However 
with the new (P or b, Mwd) correlation presented in this paper 
we obtain a larger maximum mass (i = 90°) of this pulsar: 
M$3 ax = 1.39 M© or 1.64 M© for a donor star of Pop.I or 
Pop. II chemical composition, respectively. This result brings 
the mass of PSR J20 19+2425 inside the interval of typical esti- 
mated values of Mns- 

6.4. Mns-' dependence on the propeller effect and accretion 
disk instabilities 

It is still an open question whether or not a significant amount 
of mass can be ejected from an accretion disk as a result of the 
effects of disk instabilities (Pringle 1981; van Paradijs 1996). 
However, there is clear evidence from observations of Be/X- 
ray transients that a strong braking torque acts on these neutron 
stars which spin near their equilibrium periods. The hindering 
of accretion onto these neutron stars is thought to be caused by 
their strong rotating magnetic fields which eject the incoming 
material via centrifugal acceleration - the so-called propeller 
effect (Illarionov & Sunyaev 1985). 

For a given observed BMSP we know P or b and using 
eqs (20), (21) we can find Mwd for an adopted chemical com- 
position of the donor star. Hence we are also able to calculate 
the maximum gravitational mass of the pulsar, M™g x (which 
is found for i = 90°, cf. eq. 23) since we know the mass 
function, / from observations. This semi-observational con- 
straint on M^g x can then be compared with our calculations 
of M^j S (cf. Section 5.7). The interesting cases are those where 
M^j g > M™g x (after correcting M NS for the mass deficit). 
These systems must therefore have lost matter (AMd p ^ 0), 



from the accretion disk or as a result of the propeller effect, in 
addition to what is ejected when |M 2 | > MEdd- These bina- 
ries are plotted in Fig. 4b assuming an 'intermediate' Pop. I+II 
chemical composition for the progenitor of the white dwarf. We 
notice that in a some cases we must require AMd p ~ 0.50 M©, 
or even more for M 2 > 1.0 M Q , in order to get M N s below 
the maximum limit (M^g x ) indicated by the plotted arrow. We 
therefore conclude that mass ejection, in addition to what is 
caused by super-Eddington mass-transfer rates, is very impor- 
tant in LMXBs. Whether or not this conclusion is equally valid 
for super- and sub-Eddington accreting systems is difficult to 
answer since systems which evolve through an X-ray phase 
with super-Eddington mass-transfer rates lose a large amount 
of matter from the system anyway and therefore naturally end 
up with small values of M NS . 

6.5. Kaon condensation and the maximum mass ofNS 

It has recently been demonstrated (Brown & Bethe 1994; Bethe 
& Brown 1995) that the introduction of kaon condensation suf- 
ficiently softens the equation-of-state of dense matter, so that 
NS with masses more than <~ 1.56 M© will not be stable and 
collapse into a black hole. If this scenario is correct, then we 
expect a substantial fraction of LMXBs to evolve into black 
hole binaries - unless AMd p is comparable to the difference 
between M 2 and Mwd as indicated above. However, it has re- 
cently been reported by Barziv et al. (1999) that the HMXB 
Vela X-l has a minimum value for the mass of the neutron star 
of Mns > 1-68 M© at the 99 % confidence level. It is there- 
fore still uncertain at what critical mass the NS is expected to 
collapse into a black hole. 

6.6. PSR J1603-7202 

The maximum allowed value of the pulsar mass in this sys- 
tem is extremely low compared to other BMSP systems with 
He- WD companions. We find M^ s ax = 0.96 - 1.11 M© de- 
pending on the chemical abundances of the white dwarf pro- 
genitor. It is therefore quite suggestive that this system did not 
evolve like the other BMSPs with a He-WD companion. Fur- 
thermore (as noted by Lorimer et al. 1996), it has a relatively 
slow spin period of P sp i n = 14.8 ms and P or b = 6.3 days. 
Also its location in the (P,P) diagram is atypical for a BMSP 
with a He-WD (Arzoumanian, Cordes & Wasserman 1999). All 
these characteristica are in common with BMSPs which possi- 
bly evolved through a CE evolution (van den Heuvel 1994b; 
Camilo 1996). We conclude therefore, that this system evolved 
through a phase with critical unstable mass-transfer (like in a 
CE) and hence most likely hosts a CO- WD companion rather 
than a He-WD companion. The latter depends on whether or 
not helium core burning was ignited, and thus on the value of 
P^ rh and M 2 . Spectroscopic observations should answer this 
question. 
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7. Conclusions 

- We have adapted a numerical computer code, based on 
Eggleton's code for stellar evolution, in order to carefully 
study the details of mass-transfer in LMXB systems. We 
have included, for the first time to our knowledge, other 
tidal spin-orbit couplings than magnetic braking and also 
considered wind mass-loss during the red giant stage of the 
donor star. 

- We have re-calculated the (P orrj , M^d) correlation for bi- 
nary radio pulsar systems using new input physics of stellar 
evolution in combination with detailed binary interactions. 
We find a correlation which yields a larger value of M^d 
for a given value of P orD compared to previous work. 

- Comparison between observations of BMSPs and our cal- 
culated post-accretion Mns suggests that a large amount 
of matter is lost from the LMXBs; probably as a result 
of either accretion disk instabilities or the propeller effect. 
Hence it is doubtful whether or not observations will reveal 
an (P or b , -Mns) anti-correlation which would otherwise be 
expected from our calculations. 

- The mass-transfer rate is a strongly increasing function of 
initial orbital period and donor star mass. For relatively 
close systems with light donors (P^, MS < 10 days and 
M 2 < 1.3 Mq) the mass-transfer rate is sub-Eddington, 
whereas it can be highly super-Eddington by a factor of 
~ 10 4 for wide systems with relatively heavy donor stars 
(1.6 ~ 2.0 M ), as a result of their deep convective en- 
velopes. Binaries with (sub)giant donor stars with mass in 
excess of <~ 2.0 M Q are unstable to dynamical timescale 
mass loss. Such systems will evolve through a common en- 
velope evolution leading to a short (< 10 days) orbital pe- 
riod BMSP with a heavy CO/O-Ne-Mg white dwarf com- 
panion. Binaries with unevolved heavy (> 2M Q ) donor 
stars might be dynamically stable against a CE, but also end 
up with a relatively short P orb and a CO/O-Ne-Mg WD. 

- Based on our calculations, we present new evidence that 
PSR J 1 603-7202 did not evolve through a phase with stable 
mass transfer and that it is most likely to have a CO white 
dwarf companion. 

- The pulsar mass of PSR J20 19+2425 now fits within the 
standard range of measured values for Mns, given our new 
(Porb, M WD ) correlation. 
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Appendix A: Tidal torque and dissipation rate 

We estimate the tidal torque due to the interaction between the tidally 
induced flow and the convective motions in the stellar envelope by 
means of the simple mixing-length model for turbulent viscosity v — 
aHpVc, where the mixing-length parameter a is adopted to be 2 or 
3, Hp is the local pressure scaleheight, and V c the local character- 
istic convective velocity. The rate of tidal energy dissipation can be 
expressed as (Terquem et al. 1998): 



AE 192tt 2 , 2 



dr (Al) 



where the integration is over the convective envelope and £1 is the 
orbital angular velocity, i.e. we neglect effects of stellar rotation. The 
radial and horizontal tidal displacements are approximated here by the 
values for the adiabatic equilibrium tide: 

Cr = /r2p (S 1 (A2) 

= 1 d(r 2 6Q 
6r dr 

where for the dominant quadrupole tide (I = m = 2) / = ~^Jf 2 . 
The locally dissipated tidal energy is taken into account as an extra en- 
ergy source in the standard energy balance equation of the star, while 
the corresponding tidal torque follows as: 

The thus calculated tidal angular momentum exchange dj = Ydt be- 
tween the donor star and the orbit during an evolutionary timestep dt 
is taken into account in the angular momentum balance of the system. 
If the so calculated angular momentum exchange is larger than the 
amount required to keep the donor star synchronous with the orbital 
motion of the compact star we adopt a smaller tidal angular momen- 
tum exchange (and corresponding tidal dissipation rate in the donor 
star) that keeps the donor star exactly synchronous. 
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Table Al. LMXB systems calculated for the work presented in this paper. Pop. I and Pop. II chemical compositions correspond 
to X=0.70, Z=0.02 and X=0.75, Z=0.001, respectively, a is the mixing-length parameter. P or b is in units of days and the masses 
are in units of M©. 
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10.0 


8.85 


113.5 


0.329 


2.14 


74.0 


15.0 


13.6 


162.9 


0.346 


1.85 


40.0 


25.0 


22.9 


253.1 


0.370 


1.65 


22.6 


40.0 


37.2 


369.6 


0.393 


1.54 


16.4 


60.0 


56.4 


500.5 


0.416 


1.47 


11.5 


100 


95.4 


714.6 


0.448 


1.41 


7.6 


150 


146 


933.5 


0.478* 


1.38 


5.4 


200 


197 


1108.8 


0.502* 


1.36 


4.4 



M 2 


= 1.6 


a = 


2.0 


Pop.I 




pZAMS 
-Torb 


nRLO 

-^orb 


pi 
1 orb 




Mks 


tx 


1.5 


1.09 


0.08 


0.147 


2.75 




1.8 


1.15 


1.20 


0.189 


2.71 


2350 


2.0 


1.20 


6.05 


0.219 


2.68 


780 


2.1 


1.27 


11.2 


0.246 


2.65 


450 


2.2 


1.30 


19.6 


0.263 


2.63 


285 


2.6 


1.68 


33.6 


0.283 


2.13 


136 


3.0 


2.00 


41.8 


0.291 


2.05 


123 


4.0 


2.90 


58.0 


0.303 


2.06 


107 


6.0 


5.0 


87.1 


0.320 


2.04 


78.0 


8.0 


6.8 


109.6 


0.330 


2.00 


61.0 


10.0 


8.7 


133.4 


0.339 


1.85 


42.0 


15.0 


13.3 


194.7 


0.358 


1.62 


23.2 


25.0 


22.3 


309.2 


0.384 


1.50 


14.5 


40.0 


36.0 


451.5 


0.410 


1.44 


9.3 


60.0 


54.5 


608.0 


0.435 


1.40 


7.5 


100 


92.3 


828.9 


0.466* 


1.37 


5.1 


150 


141 


1035.1 


0.494* 


1.36 


4.1 


200 


190 


1192.8 


0.531* 


1.35 


3.2 



* This is the mass of M 2 when the ~ 0.47 Mq He-core ignites (flash). 
P z r b MS is the initial orbital period of the NS and the unevolved companion. 
P ^b° is the orbital period at onset of RLO. 



Po rb is the final orbital period of the LMXB 

M theo 



or the (initial) period of the BMSP. 
WD is our calculated mass of the helium WD. 
Mn S is the final mass of the NS if AM dp = 0. 

tx is the integrated duration time (Myr) that the binary is an active X-ray source. 



M 2 


= 2.0 


Q = 


= 2.0 


Pop.I 


pZAMS 
^orb 


nRLO 

-^orb 


P f 

1 orb 


A/fthoo 
%D 




tx 


1.2 


1.19 


0.05 


0.117 


2.84 




1.3 


1.29 


2.82 


0.211 


2.75 


1820 


1.4 


1.39 


5.82 


0.232 


2.75 


1070 


1.5 


1.49 


9.56 


0.247 


2.75 


650 


1.6 


1.59 


14.9 


0.259 


2.76 


410 


1.7 


1.69 


22.7 


0.269 


2.76 


218 


1.75 


1.74 


8.1 


0.278 


2.74 


160 


1.77 


1.76 


0.5 


0.281 


2.72 


124 


1.79 


1.78 


4.1 


0.289 


2.55 


100 


1.8 


1.79 


38.6 


0.297 


1.83 


49.0 


1.9 


1.89 


39.7 


0.301 


1.82 


37.8 


2.0 


2.00 


41.0 


0.305 


1.81 


37.0 


2.2 


2.20 


44.6 


0.310 


1.79 


37.0 


2.4 


2.40 


48.2 


0.313 


1.76 


36.0 


2.6 


2.60 


51.6 


0.316 


1.74 


34.2 


2.8 


2.78 


54.7 


0.318 


1.71 


32.8 


3.0 


2.89 


56.2 


0.320 


1.64 


28.5 


3.2 


2.76 


56.4 


0.320 


1.57 


24.0 


3.4 


2.94 


59.7 


0.322 


1.57 


24.0 


3.6 


3.15 


63.3 


0.324 


1.57 


23.8 


3.8 


3.37 


66.7 


0.326 


1.57 


24.1 


4.0 


3.58 


69.8 


0.327 


1.57 


24.2 


4.2 


3.79 


73.0 


0.328 


1.57 


24.0 


>4.2 
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Table A2. Observed binary pulsars (NS+WD) in the Galactic disk. Masses are in units of M . 



PSR 


P orb (days) 


/(M Q ) 


M#5 


» ^-thco 
"WD 


i theo (deg.) 




M NS 


P sp in (ms) 


Class 


R0820+02 


1232 


0.003 


0.231 


0.503 


26.1 


6.02 


1.34 


865 


A* 


T1 803-771? 

J 1 OUJ — £u 1 L Z. 


407 


0.0013 


0.170 


0.423 


22.7 


7.20 


1.48 


334 


A 




175 


00"S8 


0.295 


0.373 


44.9 


2.61 


1.69 


3.16 


A 




147 


00078 


0.142 


0.363 


21.7 


7.48 


1.82 


4.62 


A 


B1 953+29 


117 


0.0024 


0.213 


0.352 


33.6 


3.90 


1.86 


6.13 


A 


T???Q+?|S41 


93.0 


00084 


0.146 


0.340 


23.7 


6.51 


1.90 


2.98 


A 


T20 19+2425 


76.5 


0.0107 


0.373 


0.331 


73.5 


1.51 


1.93 


3.93 


A 


TM'S'S-iiio 


76.2 


0063 


0.304 


0.331 


53.5 


2.07 


1.93 


7.99 


A 


11711+0447 


67.8 


0.0079 


0.332 


0.326 


61.6 


1.77 


1.95 


4.57 


A 


12031+1734 


56.2 


0.0027 


0.222 


0.318 


38.9 


3.13 


1.97 


5.94 


A 


R 1 8^5+09 

1 J 1 O JJTU7 


12.33 


00^7 


0.291 


0.262 


71.5 


1.54 


2.05 


5.36 


AB 


T 1804-27 17 


11.13 


00115 


0.241 


0.259 


54.0 


2.02 


2.05 


9.34 


AB 


10621+1002 


8.319 


0.0271 


0.540 


0.251 


>90 


0.51 


2.06 


28.9 


c 


11022+1001 

J ± \_7*—Z_ 1^ ± \J\J i. 


7.805 


0811 


0.872 


0.249 


>90 


0.18 


2.06 


16.5 


c 


12 145 -07*50 


6.839 


0.0242 


0.515 


0.245 


>90 


0.54 


2.07 


16.1 


c 


T2 129-5721 


6.625 


00105 

V 7 . \J\J I \ J ^ i 


0.158 


0.244 


35.4 


3.48 


2.07 


3.73 


AB 


T1 fSOl-7202 

J 1 LIU J / Z.V/Z. 


6.309 


00881 

\J.\J\JOO 1 


0.346 


0.243 


>90 


1.03 


2.07 


14.8 


Q 


10417^.7 1 5 


5.741 


001 25 


0.168 


0.240 


38.5 


3.09 


2.08 


5.76 


AB 




4 084 


001 77 


1 Q 1 


212 


A(\ 4 


2 42 


2 OQ 


7 4S 


AR 


J1911-1114 


2.717 


0.000799 


0.143 


0.222 


35.2 


3.48 


2.09 


3.63 


B 


J23 17+1439 


2.459 


0.00221 


0.206 


0.220 


54.8 


1.97 


2.10 


3.44 


B 


J02 18+4232 


2.029 


0.00204 


0.201 


0.216 


54.1 


2.00 


2.10 


2.32 


B 


B 183 1-00 


1.811 


0.000124 


0.075 


0.213 


18.8 


8.64 


2.10 


521 


B 


J0034-0534 


1.589 


0.00127 


0.169 


0.211 


45.0 


2.50 


2.11 


1.88 


B 


J06 13-0200 


1.119 


0.000972 


0.154 


0.205 


41.5 


2.78 


2.12 


3.06 


B 


B0655+64 


1.029 


0.0714 


0.814 


0.202 


>90 


0.14 


2.12 


196 


C 


J1012+5307 


0.605 


0.000580 


0.128 


0.193 


36.1 


3.33 


2.13 


5.26 


B 


B 1957+20 


0.382 


0.0000052 


0.026 


0.186 


7.3 


35.0 


2.15 


1.61 


B 


J0751+1807 


0.263 


0.000974 


0.154 


0.181 


48.2 


2.28 


2.16 


3.48 


B 


J205 1-0827 


0.099 


0.000010 


0.032 


0.168 


10.0 


21.7 


2.18 


4.51 


B 



The last column gives the classification of the BMSPs. Class A represents the wide-orbit binaries with He-WD companions. Class B 
contains the close-orbit binaries with He-WD companions. In these class B systems non-conservative angular momentum losses (J gwr , 
J m b and irradiation) were dominant in the evolution of the progenitor LMXB and P^ lb < P bi {. The subclass AB refers to systems in 
which tidal spin-orbit interactions were important but not sufficiently strong to finally prevent the orbit from widening (P^b ~ 2 — 3 
days). Class C hosts the BMSPs with heavy CO- WD companions. These systems evolved through (and survived) a phase with extreme 
mass-transfer rates and loss of orbital angular momentum (e.g. a common envelope). 



MJd is the mass of the white dwarf assuming i = 60° (the mean value of a random isotropic distribution) and Mns = 1.4 M . 
M$f? is the mass of the white dwarf obtained from our (P or b, Mwd) correlation assuming a Pop. I+II composition for the donor. 
i thco is the orbital inclination angle required for M^d = M\vfj°; i > 90° means there is no solution (the observed WD is 'too heavy'). 
M^g 51 is the maximum value for Mns obtained from the observed binary mass function, / using i — 90° and Mwd = Mwd°- 
M^g i s a rough estimate of the potential maximum post-accretion mass of the NS assuming AMd p = (see curves in Fig. 4b). 
For references to the observed values (P or b, / and P sp m) see e.g. Camilo (1995) and Lorimer et al (1996). 
* The companion of this very wide-orbit pulsar might be a CO white dwarf (cf. Section 6.2). 



